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The structures community has recognized that uncertainties in the struc- 
tural parameters as well as in the service environments need to be considered 
in evaluating structural integrity and reliability. Probabilistic structural 
analyses are formal methods to include those uncertainties. However, these 
methods are inherently computational intensive because of the large number 
of deterministic analyses required to accurately simulate the effect of 
the uncertainties on the desired structural response (such as stress, dis- 
placement, and frequency) required for structural reliability assessment. 
Moreover, modern structures are often analyzed by finite element methods. 
Probabilistic structural analyses using finite element models can be econom- 
ical when relatively coarse meshes are employed. Finite element analyses 
using coarse meshes not only raise questions regarding the convergence to the 
deterministic values, but also significantly alter the true probabilistic 
distributions of the structural responses. 

It is important, therefore, to evaluate the influence of mesh coarse- 
ness on the accuracy of structural reliability. Several alternatives were 
recently examined at Lewis. The objectives of this presentation are to 
briefly describe these alternatives and to demonstrate their effectiveness. 
The results show that special mapping methods can be developed by using 
(1) deterministic structural responses from a fine (convergent) finite ele- 
ment mesh, (2) probabilistic distributions of structural responses from a 
coarse finite element mesh, (3) the relationship between the probabilistic 
structural responses from the coarse and fine finite element meshes, and 
(4) probabilistic mapping. The structural responses from different finite 
element meshes are highly correlated. Using this correlation together with 
the probabilistic potential energy variation principle (ref. 1) that deter- 
mines the mean and standard deviation of a structural response for a given 
finite element mesh, one can obtain the linear relationships shown in 
equation (1): 
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where X f and X are the random responses from fine and coarse meshes, 
respectively, and A and B are constants. Three relationships are derived 
for three different mapping methods. Shift mapping considers only the first 
moment correction: A in equation (1) is equal to 1, and B is the 
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difference between the deterministic responses from fine and coarse meshes. 
Shift mapping states that the shapes of the probabilistic distribution func- 
tions of structural responses are the same for different meshes, but their 
mean values vary. Ratio mapping is from a second moment correction: A in 

equation (1) is the ratio of fine and coarse mesh deterministic responses, and 
B is equal to zero. Ratio mapping represents not only the fact that the mean 
values are different but also that the scatters around the mean can be either 
wider or narrower. Therefore, only the reduced response variables are con- 
sidered to be the same. Mixed mapping is the average of shift and ratio 
mapping. Once this relationship is developed, probabilistic mapping can 
be applied. 



Four examples were studied to verify this methodology. The computer code 
NESSUS (Numerical Evaluation of Stochastic Structures Under Stress) (ref. 2) 
was used to perform the probabilistic structural analyses. In the first exam- 
ple, a cantilever plate subjected to lateral pressure was analyzed. Plate 
thickness and the uniform pressure were considered to be random variables. 

In the second example, a buckling analysis of a simply supported composite 
plate was performed. The random variables were the coefficients of the stiff- 
ness matrix for the stress resultants/generalized strains relations. The 
probabilistic distributions of those coefficients were computed by the com- 
puter code PICAN (Probabilistic Integrated Composite Analyzer). In the third 
example, a cantilever plate subjected to thermal and mechanical loads was 
analyzed. Three random fields (uncertainties) were considered - thickness, 
modulus, and temperature. Each field consisted of correlated nodal random 
variables, and the loads at the free edge were also considered to be random. 
The structural parameters, such as modulus and strength, deteriorated under 
the aggressive service environments. These effects were characterized by the 
Multi-Factor Interaction Model (ref. 3). In this example, even with a poor 
mesh (16 percent error in the deterministic response), the probabilistic dis- 
tribution of the response and the structural reliability using ratio mapping 
compare very well with those using fine finite element mesh. In addition, the 
ratio of the computational time using coarse and fine mesh was about 1:20. In 
the final example, a tapered cantilever plate with variable thickness was 
studied. The uncertainties considered were the same as in the previous exam- 
ple. In each example, only the maximum effective stresses from coarse and 
fine meshes using different mapping methods were compared. From examples 1 
and 2, we verified the equality of the ratios r. and r 2 as defined in the 
viewgraph "Comparisons Between Probabilistic Distributions from Different 
Mapping Methods, Cantilever Plate Subjected to Uniform Load." This is essen- 
tial to the derivation of the ratio mapping. We also found that (1) the shift 
mapping works well only with a good coarse mesh; (2) ratio mapping, which 
provides very accurate probabilistic distribution and the structural relia- 
bility even with a very coarse mesh, is highly recoiranended; and (3) results 
from mixed mapping always lie between those from shift and ratio mapping. 


In conclusion, mapping methods were developed to perform probabilistic 
structural analyses by using coarse finite element meshes. High accuracy 
was achieved, and computational time was saved. Therefore, the dilemma 
experienced using either coarse or fine meshes for the probabilistic struc- 
tural analyses was resolved. 
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Background 


• Probabilistic structural analysis using convergent finite 
element mesh is computationally inefficient 

- Large amount of CPU time is needed 

- Long turn-around time is expected 

• Probabilistic structural analysis using coarse finite element 
mesh may not provide accurate results 

• Accuracy can be improved by 

- Shifting the mean 

- Adjusting the scatter around the mean 
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Objective 

Develop methods to increase the efficiency and 
to improve the accuracy of probabilistic structural 
analysis using coarse finite element mesh 
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Mathematical Foundation for Various 
Mapping Methods 

• Probabilistic structure responses from different finite element meshes 
are highly correlated 

• Mean and standard deviation are determined by probabilistic potential 
energy variation principle (Liu, Mani, and Belytschko (ref. 1)) 

E(X) = X° - I (X k al) 
k = 1 

o 2 = I (X k ) 2 o£ 

k =1 

where 

X probabilistic structural response for given mesh 
X° deterministic value 

Xfc 3X/dU k where U k is k ,h independent random variable 

xg a 2 x/auj; 

standard deviation of ct>-*vsji9i 


Probabilistic Relationship for Shift Mapping 

First Moment Correction 


From 

X, - E(X f ) = X c - E(X C ) 


= X c - X° 
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Since 

I (X c ) k oJ - X (X f ) k o 2 is small, 
k=1 k=1 


x,= x c -x° + xj 


• Xf and X c are the probabilistic structural responses from fine 
and coarse finite element meshes. 
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Probabilistic Relationship for Ratio Mapping 

Second Moment Correction 


From 



Since 

i-l = r 2 
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Probabilistic Relationship for Mixed Mapping 

Combined First and Second Moment Mapping 

Letting 

Xj s ~ (shift mapping) + ^ (ratio mapping) 
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Probabilistic Response Can Be Accurately 
Predicted Using Coarse Finite Element Mesh By— 

(1) Computing the convergent deterministic response 

(2) Computing the probabilistic distribution of the response 
using coarse finite element mesh 

(3) Determining the probabilistic relationship between the 
coarse and fine finite element responses 

(4) Computing the true probabilistic distribution of the 
response using the results obtained from steps (1) to (3) 
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Probability of Damage Initiation 

P f = P(c > s) 


where 

P f probability 
o stress 
s strength 


P| = r.dlMs) dsiwxldx 


f<3> fs 


Stress, a 
(from NESSUS) -7 



r Strength, s 

1 (from generic probabilistic 
1 material property model) 



Probabilistic Stress Analysis of Cantilever Plate 

• Uncertainties (random variables): plate thickness, lateral pressure 



Geometry and Load 
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Comparisons Between Probabilistic Distributions 
From Different Mapping Methods 

Cantilever Plate Subjected to Uniform Load 


Cumulative 

distribution 

function 



50 110 170 50 110 170 

Effective stress, ksi CD-91 52199 
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Comparisons Between Structural Reliabilities 
From Different Mapping Methods 

Cantilever Plate Subjected to Uniform Load 


Probability 
of damage 
initiation 
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With ratio mapping 
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Probabilistic Buckling Analysis of a Simply 
Supported Composite Plate 

• Uncertainties (random variables): coefficients in [D] matrix 
([M x My Mxy j T = [D] [Kx Ky Kxy l"^) 
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Comparisons Between Probabilistic Distributions 
From Different Mapping Methods 

Buckling Analysis of Composite Plate 


Cumulative 

distribution 

function 



2400 3200 4000 



2400 3200 4000 



Coarse mesh 


Fine mesh 
Coarse mesh 
Without mapping 
With ratio mapping 
With mixed mapping 
With shift mapping 


Buckling load, Ib/in. 
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Probabilistic Stress Analysis of Cantilever Plate 
Subjected to Thermal and Mechanical Loads 

• Uncertainties (random fields): thickness, modulus, temperature, and 
thermal and mechanical loads 




450 


1150 


Geometry and Temperature Profile (temperatures in 

Mechanical Loads degrees Fahrenheit) 
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Deterministic Stress Analyses Using Different 
Finite Element Meshes 


Nodes 

Deterministic 

stress 

Error, 

percent 

CPU time, 
sec 

16 (4 by 4) 

75.4 

16 

0.9 

25 (5 by 5) 

81.0 

9 

1.2 

36 (6 by 6) 

83.9 

6 

1.7 

49 (7 by 7) 

85.6 

4 

2.2 

64 (8 by 8) 

86.7 

3 

2.9 

81 (9 by 9) 

87.5 

2 

3.5 

100 (10 by 10) 

88.1 

1.5 

4.5 

361 (19 by 19) 

89.4 

0 

19.1 
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Comparisons Between Probabilistic Distributions 
From Different Mapping Methods 

Cantilever Plate Subjected to Thermal and Mechanical Loads 



Effective stress, ksi 



Fine mesh 


Coarse mesh 

Fine mesh 
Coarse mesh 
Without mapping 
With ratio mapping 
With mixed mapping 
With shift mapping 
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Comparisons Between Structural Reliabilities 
From Different Mapping Methods 

Cantilever Plate Subjected to Thermal and Mechanical Loads 


Probability 
of damage 
initiation 


Fine mesh 

Coarse mesh 
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A With mixed mapping 



456 7456 7 


Log fatigue cycles 


CO'S1-S2705 


Probabilistic Stress Analysis of Tapered 
Cantilever Plate 

• Uncertainties (random fields): thickness, modulus, temperature, 
and thermal and mechanical loads 



Geometry and Temperature Lateral Pressure 

CD-91-52206 


156 






Comparisons Between Probabilistic Distributions 
From Different Mapping Methods 

Tapered Cantilever Plate Subjected to Thermal and Mechanical Loads 


Cumulative 

distribution 

function 
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Comparisons Between Structural Reliabilities 
From Different Mapping Methods 

Tapered Cantilever Plate Subjected to Thermal and 
Mechanical Loads 



Fine mesh 

Coarse mesh 

O Without mapping 

□ With ratio mapping 

A With mixed mapping 
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Concluding Remarks 

Mapping methods have been developed — 

(1) To improve the accuracy of structural reliability 
using coarse finite element meshes 

(2) To save computational and turn-around time 

(3) To evaluate reliability of large structural systems 
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